Open-source Implementation of an Algorithm for Photopeaks Search and Analysis in Gamma-ray Spectrometry with Semiconductor Detectors

Radioactivity quantification of gamma-ray emitter radionuclides in samples measured by HPGe gamma spectrometers relies on the analysis of the photopeaks present in the spectra, especially on the accurate determination of their net areas. This paper presents a methodology and an algorithm description for the peak search and analysis in order to obtain the relevant peaks parameters and their uncertainties. The procedure is performed on a three step approach: a) a preliminary search is done by using the second-difference method; b) experimental peaks widths are assessed in order to obtain a width vs. channel relationship and to define regions with single or overlapping peaks; c) a non-linear fit is applied to each region of the spectrum with candidate peaks. The final target function is in the form G(x) = B(x) + F(x), where B(x) is the baseline composed by a sum of a weighed left-side BL(x) and right-side BR(x) base-line quadratic functions and the photopeaks term F(x) is a sum of Gaussian functions. The computational implementation is released entirely in open-source license. The code was developed in C++ language and the interface was developed with Qt GUI software toolkit. GNU scientific library, GSL, was employed to perform linear and non-linear fitting procedures as needed. Spectra previously generated at our laboratories were analyzed with the presented methodology and with the commercial software package WinnerGamma. Results obtained are consistent with those obtained with the aforementioned package, suggesting that it could be safely used in general-purpose gamma-ray spectrometry.